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Abstract: Determining the distributions and sources of heavy metals in soils and assessing ecological risks 
are fundamental tasks in the control and management of pollution in mining areas. In this study, we 
selected 244 sampling sites around a typical lead (Pb) and zinc (Zn) mining area in eastern Inner Mongolia 
Autonomous Region of China and measured the content of six heavy metals, including cuprum (Cu), Zn, 
Pb, arsenic (As), cadmium (Cd), and chromium (Cr). The ecological risk of heavy metals was 
comprehensively evaluated using the Geo-accumulation index, Nemerow general pollution index, and 
potential ecological risk index. The heavy metals were traced using correlation analysis and principal 
component analysis. The results showed that the highest content of heavy metals was found in 0—5 cm 
soil layer in the study area. The average content of Zn, As, Pb, Cu, Cr, and Cd was 670, 424, 235, 162, 94, 
and 4 mg/kg, respectively, all exceeding the risk screening value of agricultural soil in China. The areas 
with high content of soil heavy metals were mainly distributed near the tailings pond. The study area was 
affected by a combination of multiple heavy metals, with Cd and As reaching severe pollution levels. The 
three pathways of exposure for carcinogenic and noncarcinogenic risks were ranked as inhalation>oral 
ingestion>dermal absorption. The heavy metals in the study area posed certain hazards to human health. 
Specifically, oral ingestion of these heavy metals carried carcinogenic risks for both children and adults, as 
well as noncarcinogenic risks for children. There were differences in the sources of different heavy metals. 
The tailings pond had a large impact on the accumulation of Cd, Zn, and Pb. The source of Cr was the 
soil parent material, the source of As was mainly the soil matrix, and the source of Cu was mainly the 
nearby Cu ore. The purpose of this study is to more accurately understand the extent, scope, and source 
of heavy metals pollution near a typical mining area, providing effective help to solve the problem of 
heavy metals pollution. 
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1 Introduction 


The production of smelted lead (Pb), cuprum (Cu), and zinc (Zn) in China accounted for 50.5% of 
the global total production in 2018 (Shi et al., 2019). Chromium (Cr), Zn, Pb, arsenic (As), Cu, 
cadmium (Cd), and nickel (Ni) are listed as priority contaminants by the United States 
Environmental Protection Agency (USEPA) due to their low degradability, bioaccumulation, and 
toxicity (Muhammad, 2022). Among heavy metals, trace amounts of Cu, ferrum (Fe), and Zn are 
essential for human health. However, high concentrations of heavy metals can be toxic to humans 
and other organisms. Heavy metals, including Cd, Cr, Pb, and As, are toxic even in trace amounts 
and cause a variety of health problems. For example, Cd can cause high blood pressure, 
headaches, and osteoporosis, Cr poisoning can lead to liver and kidney dysfunction, and Pb can 
cause memory loss and anemia (Muhammad et al., 2021; Sehrish et al., 2021; Din et al., 2022). 
Mining activities are considered to be one of the main causes of heavy metals pollution in soil 
(Liu et al., 2020). The tailings pond is exposed to the natural environment over a long time and 
poses a serious threat to the entire ecosystem (Du et al., 2020; Wei et al., 2020). Heavy metals are 
not easily degradable and can accumulate in human body through the food chain and seriously 
endanger human health (Wang et al., 2016; Sun et al., 2019). 

Targeted remediation strategies can be developed by conducting ecological risk assessments 
and tracing heavy metals. Many researchers have used indices or models, such as the enrichment 
factor index, Nemerow general pollution index (Py), potential ecological risk index (RJ), and 
human health risk assessment, to simulate soil heavy metals pollution and to study the effects of 
soil physical-chemical properties on the risk of soil heavy metals pollution (Liu et al., 2022). 
Cheng et al. (2018) used the Py and RI to investigate soil heavy metals pollution in a Zn and Pb 
mining region in Yunnan Province of China, and found that 95% of the local soil is heavily 
contaminated with As. A study conducted by Chun et al. (2021) using a semi-variance-based 
model showed that the heavy metals content gradually decreased with increasing distance from 
the tailings pond in the Baiyinhua mining area of Inner Mongolia Autonomous Region, China. 
Yang et al. (2022) found that soil pH, organic matter, and cation exchange are important factors 
driving the transformation of soil heavy metals forms. 

The hazard of heavy metals depends on the total amount, distribution, chemical form, and 
transport and transformation in the environment as well as biological toxicity. Based on the 
accumulation patterns of heavy metals, ecological risk assessment is conducted to explore the 
impact of heavy metals toxicity on human health (Li et al., 2003). Zhou et al. (2015) and Sun et 
al. (2017) evaluated heavy metals pollution in mining-concentrated areas and the western mining 
areas of Inner Mongolia, respectively; they found that nonferrous metal exploration area is 
dominated by slight ecological risk, but the total carcinogenic risk exceeds the acceptable range 
for humans. In addition, insufficient knowledge of the mobility and biological effectiveness of 
heavy metals also affects the risk evaluation of the surrounding ecological environment (Dong et 
al., 2019). 

Inner Mongolia has huge mineral reserves, which is of great significance to the arrangement of 
mineral resources in China (Yang et al., 2022). Due to the fact that most of the mineral resources 
in Inner Mongolia are located in grasslands ecosystem with a weak resistance to disturbance, 
mining activities can lead to a great threat to the ecosystem. Ore mining, raw material processing 
and transportation, and tailings accumulation can cause significant soil contamination (Wang et 
al., 2018). Gao et al. (2017) found that the content of Pb, Ni, and Zn in the soil of a gold (Au) 
mining area in Xilinhot City, Inner Mongolia, is much higher than local background value. The 
result of pollution risk assessment revealed that most of the soil around a nonferrous metal mining 
area in Chifeng City, Inner Mongolia, is at moderate or heavy pollution levels, with large risk 
indices for As and Cd (Hu et al., 2018). 

An extreme lack of soil contamination assessment work in the tailings pond has affected 
local ecological remediation measures. In this study, we analyzed the pollution of Cu, Zn, Pb, 
As, Cr, and Cd at 0-20 cm soil depth in a typical Pb and Zn mining area. The spatial 
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distribution characteristics of heavy metals in the mining area were analyzed, and the pollution 
assessment and source analysis of heavy metals were conducted. This study will contribute to a 
more accurate understanding of the level and extent of heavy metals pollution, identify the 
source of heavy metals, and develop precise, economical, and efficient remediation measures. 


2 Materials and methods 


2.1 Study area 


The study area (43°20'—-43°45'N, 117°54'—118°21'E) is located in a typical Pb and Zn mining area 
of Linxi County in Chifeng City, Inner Mongolia Autonomous Region of China (Fig. 1). The area 
is topographically mountainous with an average elevation of approximately 747 m a.s.l., and the 
terrain is generally high in the northeast and low in the southwest. The climate type is semiarid 
continental monsoon with low precipitation of approximately 330-360 mm. The dominant wind 
direction is northwest, and the annual average temperature is 5.8°C. The mining area is 
approximately 1.50 km?, and the tailings pond covers an area of 0.15 km”. 


Legend Elevation (m) 
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Vy Conservatory 
Uj), Tailings pond 
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Fig. 1 Location of the study area (a) and spatial distribution of sampling sites (b) 
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2.2 Sample collection and analysis 


According to field research, topography, and wind direction, we selected the sampling sites within 
4.00 km? around the tailings pond, which are mainly distributed in farmland and woodland. We 
used Global Position System (GPS) to precisely determine the coordinates of each sampling site 
during soil sample collection. With the tailings pond as the research center and the satellite 
images as the basis, we set a 50 mx50 m grid in a 2 kmx2 km area around the tailings pond for 
sampling. Samples were collected from the top layer (0-5, 5—10, and 10—20 cm) of agricultural 
soil, and five soil samples were collected at each sampling site along the diagonal for mixing. 
Approximately 1-kg soil samples were collected using the quadrat method and brought back to 
the laboratory. A total of 732 soil samples were collected, with each layer consisting of 244 
samples. This study was conducted from October 2021 to March 2022. Soil samples were 
collected from 1 October to 8 October in 2021, and the experiment was conducted from October 
2021 to March 2022. 
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The collected soil samples were sieved to remove foreign matter such as weeds and twigs. 
After natural drying, the soil samples were ground and passed through a 100 mesh nylon sieve. 
The content of heavy metals (including Pb, Zn, Cu, As, Cd, and Cr) in soil samples was 
determined by microwave digestion with the HNO3-H2SO4-HC1O4-HF method and inductively 
coupled plasma mass spectrometry (ICAP RQplus ICP-MS, Thermo Scientific, Beijing, China). 
Three parallel samples were measured for each soil sample, and the recoveries of all elements 
after sample digestion ranged from 80% to 137%. Electrical conductivity (EC) was determined 
using the conductivity method, and pH was measured using a pH meter (STARTER 2100, 
OHAUS, Paramus, United States of America) with a ratio of soil to water of 1:5. 


2.3 Assessment methods 


The Geo-accumulation index (/geo), Py, and RI were used to evaluate the pollution of soil heavy 
metals. 

2.3.1 Geo-accumulation index (geo) 

The Jeo is used to qualitatively assess the pollution level of individual heavy metals (Muller, 
1969). It accounts for the influence of natural factors on background values, identifies the impact 
of anthropogenic activities on environment (Fei et al., 2019), and directly evaluates the pollution 
level of heavy metals through the change in actual heavy metals measurements relative to the 
environmental baseline values (Loska et al., 2004). Soil environmental quality—risk control 
standard for soil contamination of agricultural land (Ministry of Ecology and Environment of the 
People's Republic of China, 2018) was used as the baseline value for the soil pollution evaluation. 
Igeo is calculated as follows: 


C 
Lœ = lo — 1 
BET SR 2 


where C, is the measured content of heavy metals (mg/kg); and B, is the background value of 
heavy metals content in soil (mg/kg), referring to soil environmental quality-risk control standard 
for soil contamination of agricultural land quality (Ministry of Ecology and Environment of the 
People's Republic of China, 2018). The values of Jgeo are classified into seven levels, including no 
pollution (geo<0), light pollution (O</geo<1), partial moderate pollution (1</geo<2), moderate 
pollution (2</geo<3), partial heavy pollution (3<J/ge0<4), heavy pollution (4</se0<5), and severe 
pollution (Zgeo>5). 

2.3.2 Nemerow general pollution index (Py) 

The Py is an evaluation method based on a single-factor index (Liu et al., 2019). It not only 
reflects the degree of regional pollution by heavy metals but also can be used to evaluate the soil 
quality or composite pollution of whole region (Li et al., 2021). The calculation formulas of Py 
are as follows: 


P=C,/B,, 2 


2 2 
Piave) + Pimax) 
< a (3) 


where P; is single-factor index; and Piave) and Pimax) are the average and maximum values of Pi, 
respectively. We classified the Py into five levels: no pollution (P)<1), light pollution (1<P)<2), 
moderate pollution (2<Py<3), heavy pollution (3<Py<5), and severe pollution (Py=5). 

2.3.3 Potential ecological risk index (RZ) 

The RZ is characterized by its ability to reflect biological effectiveness and spatial variation 
(Hakanson, 1980), which is a commonly method for heavy metals risk assessment (Gao et al., 
2014). The calculation of RI is as follows: 


E! =T'xP, (4) 
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where ŒE; is the potential ecological risk factor of heavy metal i; and T. is the toxicity response 
factor of heavy metal i (the T. for Cu, Zn, Pb, As, Cd, and Cr is 5, 1, 5, 10, 30, and 2, 
respectively). According to the calculation results of RJ, we also classified the pollution degree of 
heavy metals, as shown in Table 1. 


Table 1 Classification of potential ecological risk 


E, Potential ecological risk level RI Potential ecological risk level 
<40 Low ecological risk <150 Low ecological risk 
40-80 Moderate ecological risk 150-300 Moderate ecological risk 
80-160 Considerable ecological risk 300-600 Considerable ecological risk 
160-320 Very high ecological risk >600 Very high ecological risk 
>320 Serious ecological risk 


Note: E, potential ecological risk factor of a single heavy metal; RZ, potential ecological risk index. 


2.3.4 Human health risk assessment 

In this study, we consulted the Exposure Factors Handbook (U.S. Environmental Protection 
Agency, 2011) to assess carcinogenic and noncarcinogenic risks of heavy metals (Young et al., 
2022). There are three ways humans are exposed to heavy metals: oral ingestion, inhalation, and 
dermal absorption. The equations for the average daily exposure, noncarcinogenic risk, and 
carcinogenic risk are provided as follows: 


C, x [Ring x EF x ED x10 


ADD. „= , 6 
i BWxAT (6) 
‘Lx SAx ABS x EF x EDx10° 
P Wee ae x SAx Sx x x10 (7) 
BW x AT 
ADD gs C, x IRn X ED x EF , (8) 
BW x PEF x AT 


where ADDing is the average daily intake of heavy metals by oral ingestion (mg/(kg-d)); JRing is 
the ingestion rate by oral ingestion (mg/d); EF is the frequency of exposure (d/a); ED is the 
exposure duration (a); BW is the body weight (kg); AT is the average exposure time (d); ADDaer is 
the average daily intake of heavy metals by dermal absorption (mg/(kg-d)); SZ is the adherence 
factor (mg/(cm?.d)); SA is the skin exposure area (cm’); ABS is the skin exposure factor; ADDinh is 
the average daily intake of heavy metals by inhalation (mg/(kg-d)); /Rinh is the ingestion rate by 
inhalation (m?/d); and PEF is the peak expiratory flow (m°/kg). 


ADD 
HI=)\HO=>° RD’ (9) 
TCR=)\CR=)\ ADDx SF , (10) 


where HO and CR are the noncarcinogenic and carcinogenic risks for a single heavy metal, 
respectively; HI and TCR are the noncarcinogenic and carcinogenic risks for multiple heavy metals, 
respectively; ADD is the daily average exposure dose for noncarcinogenic risks (mg/(kg-d)); RD is 
the reference dose (mg/(kg-d)); and SF is the safety factor. According to the Exposure Factors 
Handbook (U.S. Environmental Protection Agency, 2011), we concluded that when CR<10°, heavy 
metals pose little or no hazard to humans; when 10°<CR<10~, heavy metals are in the acceptable 
range to humans; and when CR>10~, heavy metals pose an unacceptable risk to humans. 


2.4 Statistical methods 


We used SPSS 22.0 software to conduct descriptive statistics of soil heavy metals content and soil 
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physical-chemical properties. Correlation analysis of the heavy metals content with soil 
physical-chemical properties was performed using SPSS 22.0 software. OriginPro software was 
used to analyze the heavy metals content and soil physical-chemical properties at different soil 
depths. The spatial distribution of heavy metals risk was mapped using ArcGIS 10.2 software and 
the ordinary kriging method. Correlation analysis, principal component analysis (PCA), and 
cluster analysis were used to identify the sources of heavy metals in soils. 


3 Results 


3.1 Descriptive statistics of soil heavy metals 


The results of heavy metals content in the study area are shown in Figure 2. The average content 
of Zn, As, Pb, Cu, Cr, and Cd at 0-5 cm soil depth was 671, 425, 235, 163, 95, and 4 mg/kg, 
respectively, which were higher than the risk screening values for soil contamination of 
agricultural land required by soil environmental quality—risk control standard for soil 
contamination of agricultural land (Ministry of Ecology and Environment of the People's 
Republic of China, 2018). These results also exceeded the local (Chifeng City) soil heavy metals 
background value. The exceeding standard rate of As, Cd, Zn, Cu, Pb, and Cr was 78%, 73%, 
58%, 53%, 34%, and 10%, respectively, compared with the risk screening values (Table S1). The 
mean content of Zn, As, Pb, Cu, Cr, and Cd at 5-10 cm soil depth was 660, 354, 200, 162, 88, and 
3 mg/kg, respectively, which exceeded the background values by 18, 45, 11, 11, 2, and 64 times, 
respectively. The content of As, Cd, Zn, Cu, Pb, and Cr at 5-10 cm soil depth exceeded the risk 
screening values by 73%, 66%, 56%, 49%, 29%, and 6%, respectively. At 10-20 cm soil depth, 
the content of As, Cd, Zn, Cu, Pb, and Cr exceeded the background values by 17%, 32%, 9%, 
10%, 2%, and 104%, respectively, and the risk screening values by 64%, 62%, 52%, 41%, 23%, 
and 5%, respectively. The maximum exceeding standard rate of As and Cd was most severe at 
different soil depths, and the exceeding standard rate decreased continuously with increasing soil 
depth. The pollution level of Cr was relatively light. 

The coefficient of variation (CV) is an important indicator of data dispersion and is graded as 
weak variability (<<10%), moderate variability (10%—100%), and strong variability (>100%). The 
CV of Cr content was 87% and showed a moderate variability (Table S1). The CV values of other 
heavy metals showed strong variability, with As showing the largest CV (220%). 
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Fig. 2 Content of cuprum (Cu; a), lead (Pb; b), zinc (Zn; c), arsenic (As; d), chromium (Cr; e), cadmium (Cd; f), 
at different soil depths. The bar chart and the width of violin chart represent the frequency of the data. The black 
dots represent the medians. The bars represent the 95% confidence interval. 
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3.2 Soil physical-chemical properties 


Soil pH and EC affect a variety of soil properties, which in turn affect the biological activity of 
soil heavy metals. Results of the statistical analysis of soil physical-chemical properties are 
shown in Figure 3. At 0-5 cm soil depth, the EC ranged from 57 to 3098 mS/cm, with a mean 
value of 415 mS/cm; the mean value of soil pH was 7.7, and 16 soil samples had a pH value lower 
than 6.5. At 5-10 cm soil depth, the EC ranged from 37 to 2841 mS/cm, with a mean value of 455 
mS/cm; the value of soil pH varied from 2.4 to 8.5, with a mean value of 7.7. At 10—20 cm soil 
depth, the mean value of EC and pH was 557 mS/cm and 7.6, respectively. Most of the soil 
samples in the study area were weakly alkaline, and with the deepening of the soil depth, pH 
gradually decreased and EC gradually increased. 
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Fig. 3 Soil pH (a) and electrical conductivity (EC; b) at different soil depths. The bar chart and the width of 
violin chart represent the frequency of the data. The black dots represent the medians. The bars represent the 95% 
confidence interval. 


3.3 Spatial distribution of the content of heavy metals in soils 


The spatial distribution of the content of heavy metals in soils of the study area was obtained by 
the ordinary kriging method (Fig. 4). The distributions of the six heavy metals (Cu, As, Pb, Zn, 
Cr, and Cd) were very similar, showing a downward trend from southwest to northeast. High Cr 
contents were concentrated in the west of the tailings pond, while high As contents were 
concentrated in the southeast of the tailings pond. The highest content of heavy metals occurred at 
0-5 cm soil depth and the lowest content of heavy metals occurred at 10-20 cm soil depth (Figs. 
S1 and S2). 


3.4 Spatial distribution of heavy metals pollution in soils 


3.4.1 Spatial distribution of Zgeo 

A map of the spatial distribution of soil heavy metals pollution was drawn using the ordinary 
kriging method (Fig. 5). Results showed that the pollution level near the tailing ponds is more 
severe. The spatial distribution characteristics of different heavy metals greatly varied. The 
pollution degree of heavy metals tended to decrease with increasing soil depth. The 0-5 cm soil 
layer was the most polluted. The highest As pollution index was found at 0—10 cm soil depth in 
the southern and eastern parts of the tailings pond, with severe pollution (/geo>5). Cd pollution 
was more serious in the southwest of the tailings pond, with geo ranging from 2 to 5. 

In the study area, 7% of the soil was severely polluted by As, and 10% of the soil was 
moderately polluted by As. Additionally, 2% of the soil was heavily polluted by Cd, while 6% of 
the soil was lightly polluted by Cd (Fig. S3). The pollution degree of Cu, Pb, and Zn around the 
tailings pond was light. The soil in the western and southern parts of the tailings pond was 
moderately polluted by Cu (2</geo<3). In the northwest and southeast of the tailings pond, soil 
was moderately polluted by Zn (2</ge0<3). Soil contaminated by Zn, Cu, and Pb accounted for 
9%, 3%, and 5% of the moderately polluted areas, respectively, and 2%, 1%, and 1% of the 
heavily polluted areas, respectively. 
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Fig. 4 Spatial distribution of the content of Pb (a), Zn (b), Cu (c), As (d), Cr (e), and Cd (f) at 0-5 cm soil depth 


3.4.2 Overall pollution levels of heavy metals based on the Nemerow general pollution index (Py) 


As, Cd, Pb, Cu, and Zn reached severe pollution levels in the study area. The pollution level of As 
was the most severe, with Py ranging from 127 to 250 (Fig. 6). The pollution level of Cd was also 
more severe, with Py greater than or equal to 41. The pollution level of Cr was light, with a Py 
value of 1 at 0-5 cm soil depth. Py was less than 1 at 5-20 cm soil depth, indicating a state of no 
pollution. In summary, the pollution level of As and Cd in the study area was severe, and the 
pollution level of Cr was light. 


3.4.3 Comprehensive pollution status of heavy metals based on the potential ecological risk 
index (RZ) 

We used the soil environmental quality—trisk control standard for soil contamination of 
agricultural land (Ministry of Ecology and Environment of the People's Republic of China, 2018) 
as the baseline value to evaluate soil pollution in the study area. The ordinary kriging method was 
used to draw the spatial distribution of RZ at different soil depths in the study area (Fig. 7). 

The values of potential ecological risk factor of a single heavy metal (E) of Cr, Zn, Cu, and Pb 
at different soil depths in the study area were lower than 40, which indicated that these heavy 
metals are at a low ecological risk level (Table S2). The Æ. value of As ranged from 95 to 160, 
indicating a considerable ecological risk. The Æ. value of Cd ranged from 123 to 209, showing a 
very high ecological risk. The percentage of Cd causing serious ecological risk to soil was greater 
than that of As (Fig. S4). Therefore, the ecological risk of Cd pollution was higher than that of As 
pollution. The RJ was 388, 273, and 291 at 0-5, 5-10, and 10-20 cm soil depths, respectively. 
The 0-5 cm soil layer was at a considerable ecological risk, and the 5—20 cm soil layer was at a 
moderate ecological risk (Table S2). 
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Fig. 5 Spatial distribution of heavy metals pollution level at 0-5, 5—10, and 10-20 cm soil depths based on the 
Geo-accumulation index (/geo). (al—a3), Pb; (b1—b3), Zn; (cl—c3), Cu; (d1—d3), As; (el—e3), Cd; (f1—-f3), Cr. 


As showing in Figure 7, the heavy metals pollution at 0-5 cm soil depth was more severe in the 
southern area of the tailings pond. There was a very high ecological risk in a small part of the 
western and northern areas of the tailings pond. The areas with very high ecological risk and 
considerable ecological risk accounted for 20% and 13% of the total area, respectively. The 
pollution degree at 5—20 cm soil depth gradually decreased, and the heavily polluted areas were 
mainly concentrated in the southern part of the study area. The area with very high ecological risk 
at 5-10 cm soil depth was 11%, which was 8% lower than the ecological risk of topsoil (0-5 cm). 
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The area with considerable ecological risk at 5—10 cm soil depth was 13%, showing no significant 
variation compared to the topsoil (0-5 cm). The areas with very high ecological risk and 
considerable ecological risk at 10-20 cm soil depth dropped to 7% and 10% of the total area, 
respectively (Fig. S5). The results of the RJ showed that the pollution level of the topsoil (0-5 cm) 
was the most serious, and the RZ showed a decreasing trend with the deepening of the soil layer. 


3.4.4 Health risk assessment 


The noncarcinogenic health risks of five heavy metals, including Cu, Zn, As, Pb, and Cd, were 
assessed in this study. The results showed that As, Zn, and Pb pose high noncarcinogenic risks to 
human health (Table 2). The ranking of heavy metals intake levels through different pathways was 
as follows: HOing>HOder>HOinn (where HQing HOder, and HOinn are noncarcinogenic risk of a 
single heavy metal through oral ingestion, dermal absorption, and inhalation, respectively). The 
noncarcinogenic risks for multiple heavy metals (H7) for children was 1.2, indicating that heavy 
metals in the study area pose a potential noncarcinogenic health risk threat to children. All other 
exposure pathways had an HI value of less than 1.0, indicating no significant risk to human 
health. 


Table 2 Assessment of reference dose (RfD) values and noncarcinogenic health risks of different exposure 
pathways of soil heavy metals in different populations 


Heavy RfDing  RfDar RYDinn HOn Des scion 

metal mg/(kgd) mg/(kg:d) mg/(ke-d) Child Adult Child Adult Child Adult 
Cu 40x10?  1.2x10? - 1.6x10  9.2x105  3,8x107  3.2x107  3.8x10° 8.6x10° 
Zn 3.0x10  6.0x102 g 6.4x104 38x104 16x10%  1.3x10%  1.5x105  3.5x105 


As 3.0x10* 1.2x104 3.5x10°  3.4x10* 2.0x104 8.3x107  7.0x107 8.2x10° 1.8x10® 
Pb 3.5x103 5.3x104 8.2x105 1.8x104 1.0x104 4.3x107  3.6x107 = 4.3x10° = 9.6x10°° 
Cd 1.0x10% 1.0x10° = 2.4x10° ~=—-2.3x10°° 1.4x10°° 5.7x10°  4.8x10° = 5.6x107'!! 1.3x107"° 


HI - - - 1.2 7.0x107! 8.2x10°  6.9x10° 2.4x10°  5.4x10° 


Note: RfDing, reference dose of a single heavy metal through oral ingestion; RfDaer, reference dose of a single heavy metal through 
dermal absorption; RfDinn, reference dose of a single heavy metal through inhalation; HQjng, noncarcinogenic risk of a single heavy 
metal through oral ingestion; HQaer, noncarcinogenic risk of a single heavy metal through dermal absorption; HQinh, noncarcinogenic 
risk of a single heavy metal through inhalation; HI, noncarcinogenic risks for multiple heavy metals. -, no value. 

The carcinogenic health risk of heavy metals was ranked as CRing>CRaer>CRinn (where CRing 
CReer, and CRinh are carcinogenic risk of a single heavy metal through oral ingestion, dermal 
absorption, and inhalation, respectively). As and Pb posed a carcinogenic risk to human health, 
while the carcinogenic risk of Ca was within an acceptable range. Dermal absorption and inhalation 
to heavy metals did not pose a carcinogenic risk to human health (Table 3). According to the 
Exposure Factors Handbook (U.S. Environmental Protection Agency, 2011) for children and adults, 
we found that oral exposure to heavy metals has an unacceptable carcinogenic risk, dermal exposure 
to heavy metals has an acceptable risk, and inhalation of heavy metals has no carcinogenic risk. 


Table 3 Assessment of carcinogenic health risks of different exposure pathways of soil heavy metals in different 
populations 


Heavy CRing CReer CRinn 

metal Child Adult Child Adult Child Adult 
As 3.4x104 2.0x10* 8.3x107 7.0x107 8.2x10° 1.8x10°8 
Pb 1.8x10~ 1.0x10~ 4.3x107 3.6x107 4.3x10° 9.6x10° 
Cd 23x10 1.4x10°° 5.7x10° 4.8x10° 5.61071! 1.3x10-! 


Note: CRing, carcinogenic risk of a single heavy metal through oral ingestion; CRaer, carcinogenic risk of a single heavy metal through 
dermal absorption; CRinn, carcinogenic risk of a single heavy metal through inhalation. 
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3.5 Multivariate statistical analysis results 


3.5.1 Correlation analysis among heavy metals 

Correlation analysis is an important basis for determining the source of heavy metals. The 
normality test was performed by SPSS software, and all the six heavy metals were non-normal 
distribution with outliers. Therefore, the Spearman correlation analysis was used in this study 
(Table 4). Cd showed highly significant positive correlations with Pb (7=0.84), Cu (7=0.58), and 
As (r=0.49); Cu showed significant positive correlations with Zn (r=0.35) and Pb (7=0.53); Pb 
showed positive correlations with As (7=0.49) and Zn (7=0.33); and Cd showed a positive 
correlation with Zn (r=0.41). There were positive correlations among Cu, Zn, Pb, Cd, and As, 
suggesting that they may have similar origins. Cr was significantly positively correlated with Pb, 
while Cr was significantly negatively correlated with Zn. The positive correlation between Cr and 
other heavy metals was not significant, indicating that the sources of Cr and other heavy metals 
are different. In addition, the pH was negatively correlated with Cu, Pb, Zn, As, and Cd and 
positively correlated with Cr. Among them, pH, EC, and Cu were significantly correlated with 
Zn; therefore, soil pH and EC in the study area may have a greater influence on the enrichment of 
Cu and Zn in soils. 


Table 4 Spearman correlation analysis of heavy metals and soil physical-chemical parameters 


Cu Zn Cr As Cd Pb EC pH 
Cu 1.000 
Zn 0.3517 1.000 
Cr 0.103 -0.355" 1.000 
As 0.196% 0.199% —0.097 1.000 
Cd 0.583” 0.411" 0.116 0.485" 1.000 
Pb 0.532" 0.331" 0.214" 0.491 0.8407 1.000 
EC 0.382" 0.366" —0.134" 0.138” 0.282" 0.229" 1.000 
pH -0.202™ -0.370" 0.226" —0.045 0.115 0.059 -0.476" 1.000 


Note: EC, electrical conductivity; *, P<0.05 level; **, P<0.01 level. 


3.5.2 Principal component analysis (PCA) 

Furthermore, in order to investigate the interrelationships among heavy metals and the possible 
sources of heavy metals, we conducted a PCA using SPSS software. The results of PCA are 
presented in Table 5. 


Table 5 Principal component analysis (PCA) of heavy metals in soils of the study area 


Initial eigenvalue Sum of squared loadings 
PC Hisseruiig a e e ii Eigenvalue ne een Cumulative variance (%) 
1 2.583 43.054 43.054 2.583 43.054 43.054 
2 1.436 23.937 66.991 1.436 23.937 66.991 
3 0.814 13.566 80.557 0.814 13.566 80.557 
4 0.575 9.580 90.137 
3 0.327 5.449 95.586 
6 0.265 4.414 100.000 


Note: PC, principal component. 


From the results of PCA (Table 5), we can see that the first three principal components of the 
six heavy metals explained 80.557% of the total variance, so the analysis of the first three 
principal components provides most information about the contents of Cu, Zn, As, Pb, Cd, and Cr 
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(Fig. 7). The PC1 (PC is principal component) explained 43.054% of the total variance, mainly 
contributed by Pb, Cd, and Zn, with rotation factor loadings of 0.87, 0.86 and 0.64, respectively. 
The PC2 explained 23.937% of the total variance, mainly contributed by Cr and As, with rotation 
factor loadings of 0.81 and 0.51, respectively. The PC3 explained 13.566% of the total variance, 
mainly contributed by Cu, with a rotation factor loading of 0.84. The results indicated that there are 
three main sources of soil heavy metals in the study area. 


PC2 (23.937%) 
Oo 
T 


Fig. 8 Principal component analysis (PCA) of heavy metals in soils of the study area. PC, principal component. 


3.5.3 Cluster analysis 

Hierarchical clustering dendrogram can directly reflect the correlation among heavy metals and 
reveal the sources of soil heavy metals. Results of the cluster analysis (the Ward method was 
chosen for the cluster analysis, and the Euclidean distance method was chosen for the 
measurement of interval) were used to group the six heavy metals into three categories (Fig. 8). 
The first group included Cd, Pb, and Zn, the second group included Cr and As, and the third 


group included Cu. 
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Cu 
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Fig.9 Cluster analysis of heavy metals in soils of the study area 


4 Discussion 


4.1 Effects of heavy metals pollution risk 
Heavy metals are found in soil during natural formation and are generally harmless to the 
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environment due to their low content (Muhammad et al., 2022). However, with the increase of the 
frequency of human activities, the content of heavy metals in the soil continues to increase, 
resulting in the content of heavy metals in the affected soil being much higher than the natural 
background value. Therefore, heavy metals pose a serious risk to ecosystems and human health 
(Li et al., 2018). 

In this study, areas with a higher degree of soil heavy metals pollution and ecological risk are 
concentrated in the northwest and southeast of the tailings pond. Therefore, heavy metals 
generally migrate from northwest to southeast. The study area features a rugged topography with 
significant variations in elevation, gradually decreasing from northwest to southeast. The 
distribution of heavy metals is more consistent with wind direction and topography, which 
indicates that wind direction and topography have an influence on heavy metals migration. This 
promotes heavy metals particulate matter entering the surrounding soil through the air. However, 
the heavy metals content decreases with increasing soil depth (Figs. 4, S1, and S2). Heavy metals 
such as Cd, As, and Zn mainly accumulate on the surface, indicating that the soils may be heavily 
polluted by exogenous sources such as atmospheric deposition and industrial and agricultural 
activities (Song et al., 2018). The heavy metals content in soils far away from the tailings pond is 
significantly lower than that around the tailings pond, and the vertical change of the heavy metals 
content is small (Zhao et al., 2021). 

Results of the correlation analysis between soil pH and heavy metals provide another concern. 
Soil pH can regulate activities related to the toxicity and transformation of heavy metals in soils 
(Xu et al., 2017). In this study, pH is negatively correlated with all heavy metals except Cr, which 
indicates that the lower the pH, the higher the heavy metals content. The effectiveness of heavy 
metals in soil solution increases with the decrease of pH (Inka et al., 2016). Therefore, the 
increased toxicity of heavy metals caused by soil acidification requires special attention. 


4.2 Analysis of the sources of heavy metals 


Human activities have an impact on the content of Cu, Zn, Cd, Pb, Cr, and As in soils. A 
combination of multivariate statistical analysis and spatial analysis was used to analyze the 
sources of heavy metals. 

Both correlation analysis and PCA suggest that Cd, Pb, and Zn may have a common source. 
Notably, the main source of Pb and Zn is the accumulation of wastewater and sludge in the 
tailings pond (Ma et al., 2015). Results of the spatial distribution of heavy metals in this study 
show that Pb and Zn are mainly concentrated near the tailings pond with characteristics of 
diffusion to surrounding soils. This result indicates that Pb and Zn are more influenced by the 
tailings pond. The regions with high Zn content are mainly located in the northwest of the tailings 
pond. The mining area processes and refines the crude ore into concentrate. This process 
generates large quantities of toxic and hazardous substances containing Zn, and the release of 
aerosols and dust can contaminate nearby soils (Liu et al., 2018). It is inferred that the main 
source of Cd is the accumulation of wastewater and slag. Cd is extremely harmful to humans and 
has a very low threshold compared to other heavy metals. Under the same conditions, Cd has 
higher migration capacity than Pb and Zn, and is more likely to contaminate deep soils and 
groundwater. As a result, the pollution level of Cd in soils of the study area is higher. 

The PC2 includes Cr and As, with most of the Cr accumulation in the study area close to the 
soil background value. Only a few soils close to the tailings pond were lightly contaminated by 
Cr. The main source of Cr in soils is the soil parent material (Liang et al., 2017; Liu et al., 2021). 
In addition, the soil in the study area generally exhibits high levels of As content. Studies have 
shown that Inner Mongolia is the region with the highest As emission potential in China (Wei et 
al., 2020). With the passage of time, the background value of As in the study area gradually 
increases. Combined with the correlation analysis, it is clear that Cr and As are homologous. 
Therefore, the main source of Cr is the soil parent material, and the study area has a high 
background value of As. 

The PC3 includes Cu. As shown in Figure 4, the areas with high Cu content are concentrated in 
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the northwest of the tailings pond. The investigation found a Cu mine approximately 1 km 
northwest of the tailings pond. The mining and smelting of Cu produces large amounts of 
elemental Cu, which can accumulate in surrounding soils through atmospheric deposition (Nezat 
et al., 2017). The dominant wind direction in the study area is northwester; therefore, the main 
source of Cu is likely the Cu mines in the northwestern part of the study area. 

Additionally, most chemical fertilizers and pesticides contain a certain amount of heavy metals. 
Even if chemical fertilizers and pesticides meet national standards, agricultural soil will be 
continuously enriched with heavy metals in the process of long-term application (Fei et al., 2019). 


5 Conclusions 


In this study, we investigated the heavy metals pollution in soils around a typical Pb and Zn mining 
area in eastern Inner Mongolia. The content and distribution of Zn, As, Pb, Cu, Cr, and Cd were 
analyzed. The ecological risk of heavy metals was assessed, and a traceability study was conducted. 
The most severely polluted heavy metals were found at 0-5 cm soil depth. The average content of 
Zn, As, Pb, Cu, Cr, and Cd was 670, 424, 235, 162, 94 and 4 mg/kg, respectively. There were 
significant differences in heavy metals contents in different regions, which were mainly influenced 
by topography and wind direction. Pollution assessment showed that the ecological risk of heavy 
metals near the tailings pond is the highest, and the impact range is the largest with some areas 
reaching severe pollution levels. The human health risk assessment showed that the three pathways 
of exposure for carcinogenic and noncarcinogenic risks are ranked as inhalation>oral 
ingestion>dermal absorption. The noncarcinogenic risks for As, Pb, and Zn were high, and in terms 
of total noncarcinogenic risk, oral ingestion posed a potential threat to children, whereas the other 
pathways did not pose a significant threat to humans. The overall carcinogenic risk was high for oral 
ingestion, within acceptable limits for dermal absorption, and did not pose a carcinogenic risk for 
inhalation. Analysis of the sources of heavy metals showed that the tailings pond contributed to the 
pollution of Pb, Zn, and Cd. The source of Cr was the soil parent material, while As was mainly due 
to high background values in the study area; and the source of Cu was mainly from nearby Cu 
mining and smelting activities. The remediation strategies in this region should focus on controlling 
the emissions of heavy metal at source. Combined with different remediation strategies, economical, 
rapid, and efficient soil remediation actions should be taken. 
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Fig. S1 Spatial distribution of the content of lead (Pb; a), zinc (Zn; b), cuprum (Cu; c), arsenic (As; d), chromium 


(Cr; e), and cadmium (Cd; f) 
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Fig. S2 Spatial distribution of the content of Pb (a), Zn (b), Cu (c), As (d), Cr (e), and Cd (f) at 10-20 cm soil depth 
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Fig. S3 Percentage of the geo-accumulation index (/geo) at 0-5 (a), 5—10 (b), and 10-20 (c) cm soil depths 


(a) 0-5 cm (b) 5-10 cm 
100- 1% 1% 100. 2% 2% 
804 
‘3 60+ D 
F E 
al E 
À A 
20+ 
0 0. 
Cd Zn Cu Pb Cr As Cd Zn Cu Pb Cr 
Heavy metal Heavy metal 
(c) 10-20 cm 
100- 1% 1% 
805 
a a Serious ecological risk 
5 60-4 sl Very high ecological risk 
o 
2 O | Considerable ecological risk 
5 404 EL] Moderate ecological risk 
g [| Low ecological risk 
20+ 
0 
As Cd Zn Cu Pb Cr 
Heavy metal 


Fig. S4 Percentage of the potential ecological risk factor of a single heavy metal (Ei) at 0-5 (a), 5-10 (b), and 
10-20 (c) cm soil depths 
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Fig. S5 Percentage of potential ecological risk index (RJ) at 0-5 (a), 5—10 (b), and 10-20 (c) cm soil depths 
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Table S1 Statistical analysis results of soil heavy metal content 


Soil depth Parameter Cu Pb Zn As Cr Cd pH EC Elevation 
(cm) (mg/kg) (mg/kg) (mg/kg) (mg/kg) (mg/kg) (mg/kg) (uS/cm) (m) 
Mean 163 235 671 425 95 4 77 415 747 
Median 105 101 366 119 64 2 80 213 748 
Minimum 3 0 18 6 9 0 2.5 57 703 
9-5 Maximum 1025 2943 5654 8840 462 87 85 3098 788 
SD 163 373 841 933 83 8 0.9 482 20 
CV (%) 100 159 125 220 87 187 117 116 3 
pais standard: -53 34 58 78 10 B 
Mean 162 200 200 354 88 3 7 455 747 
Median 86 80 358 105 59 1 7.9 235 748 
Minimum 1 1 8 2 11 0 24 37 703 
5—10 Maximum 2776 1757 6037 6458 306 20 86 2841 788 
SD 258 294 886 741 67 3 10 488 20 
CV (%) 159 147 444 209 76 125 124 107 3 
“ri standard 49 29 56 B 6 66 
Mean 138 168 620 257 91 4 76 557 747 
Median 73 70 304 61 62 1 7.9 302 748 
Minimum 0 1 0 2 10 0 27 33 703 
10-29 Maximum 1924 1497 5207 5985 2015 416 8.8 2713 788 
SD 198 257 828 656 142 30 09 546 20 
CV (%) 144 153 134 256 155 722 123 98 3 
ar standard 41 23 52 64 5 62 


Note: CV, coefficient of variation. 


Table S2 Mean value of E.and RI of soils at different depths 


Soil depth (cm) i RI 
Cu Zn Cr As Cd Pb 

0-5 10 2 1 160 209 7 388 

5-10 9 2 1 132 123 6 273 

10-20 8 2 1 95 182 4 291 


Note: E;, potential ecological risk factor of single heavy metal; R7, potential ecological risk index. 


